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We study the spin-1/2 two-leg Heisenberg ladder with four-spin ring exchanges under a 
magnetic field. We introduce an exact duality transformation which is an extension of the 
spin-chirality duality developed previously and yields a new self-dual surface in the parameter 
space. We then determine the magnetic phase diagram using the numerical approaches of the 
density-matrix renormalization-group and exact diagonalization methods. We demonstrate 
the appearance of a magnetization plateau and the Tomonaga-Luttinger liquid with dominant 
vector-chirality quasi-long-range order for a wide parameter regime of strong ring exchange. 
A "nematic" phase, in which magnons form bound pairs and the magnon-pairing correlation 
functions dominate, is also identified. 

KEYWORDS: ring exchange, spin ladder, frustration, vector chirality, nematic phase, density- 
matrix renormalization-group method 

1. Introduction 

Multiple-spin ring-exchange process - a process in which several particles with spin per- 
mute their positions in a cyclic fashion - has been found to be important in a wide vari- 
ety of materials. It has been established that ring exchanges are responsible for the mag- 
netism of solid ^He^' as well as of two-dimensional quantum solids including ^He adsorbed 
on graphite^"^^ and Wigner crystals.^' ^"'^^ Ring exchanges were also found to be relevant 
to some strongly correlated electron systems such as the two-leg spin ladder compound 
LaxCai4_xCu2404i^"^^-' and two-dimensional antiferromagnet La2Cu04.^^"^^^ It has been sug- 
gested recently that ring exchanges can be dominant in quantum wires with a shallow confining 
potential.^"' ^i) 

In contrast to two- and three-spin exchange processes, the ring exchanges involving more 
than three spins take a very complicated form in terms of spin operators including bilin- 
ear, biquadratic, and higher-order couplings in general, which compete with each other. The 
multiple-spin ring exchanges therefore contain strong frustration in themselves and can in- 
duce exotic phenomena. Thus motivated, extensive studies have been devoted to clarifying the 
eff'ects of ring exchanges on quantum spin systems in two-leg ladder, ^^"^^^ square, ^^"^^^ and 
triangular lattices. '^^"^^^ It has been revealed that the four-spin ring exchange can lead to novel 
quantum states, e.g., the scalar-chirality ordered state in the two-leg ladder system, ^^"^'^^ the 
vector-chirality^'^^ and nematic^^^ orderings in the square lattice, and the octapolar order^^^ 
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in the triangular lattice. 

In this paper, we study the S = 1/2 two-leg Heisenberg ladder system with four-spin ring 
exchanges under a magnetic field. The magnetic field lowers the symmetry of the model from 
SU(2) to U(l), and as a result, opens a possibility of novel quantum states. The model was 
studied previously by numerical exact-diagonalization method and perturbation theory and it 
was predicted that the model exhibited a plateau in the magnetization curve^"^'^^) and a field- 
induced incommensurate quasi long-range order (LRO) around the plateau. '^^^ Further, it was 
also suggested very recently that the interplay between four-spin exchanges and the magnetic 
field might give rise to a true LRO of vector chirality.^*^^ Nevertheless, these studies were 
limited to antiferromagnetic bilinear couplings and rather weak ring exchange, and therefore, 
the properties of the model for the regime of strong ring exchange and/or of ferromagnetic 
bilinear coupling remain unclear. 

The aim of the paper is to investigate the ground-state properties of the ladder system 
with four-spin ring exchanges for the whole parameter space of both the antiferromagnetic and 
ferromagnetic bilinear exchanges and a finite magnetic field. We first employ a theoretical ap- 
proach using duality transformations. The so-called spin-chirality duality transformation^^' ^'''^ 
has proven to be a powerful tool to study the ladder system with four-spin exchanges. We fur- 
ther introduce another duality transformation, which is a simple extension of the spin-chirality 
duality but gives a nontrivial duality mapping on coupling constants in the Hamiltonian as 
well as on various order parameters. In particular, applying the latter transformation to the 
ladder with the ring exchange yields a new self-dual line, which divides the parameter region 
of strong ring exchange from that of strong ferromagnetic bilinear exchange. 

Moreover, we study the system using numerical approaches of the density-matrix 
renormalization-group (DMRG) and exact diagonalization methods and determine the mag- 
netic phase diagram. We then directly observe the magnetization plateau appearing for a 
wide parameter region of strong ring-exchange coupling. We find various critical states with 
different dominant quasi-LRO: In addition to the critical state with dominant Neel-type- 
spin quasi-LRO found in usual antiferromagnetic two-leg ladder without ring exchange, the 
vector-chirality-dominant critical state appears for strong ring exchange. More interestingly, 
we show that a nematic phase, in which magnons form bound pairs and the magnon-pairing 
correlations are dominant, emerges in a finite region adjacent to the ferromagnetic phase. 

The paper is organized as follows. In Section 2, we introduce the model Hamiltonian and 
summarize its known properties obtained from previous studies. The duality transformations 
and their application to the ladder system with four-spin exchanges are discussed in Section 
3. Section 4 shows our numerical results on the magnetic phase diagram. The results on the 
magnetization plateau and the vector-chirality-dominant critical phase appearing in the region 
of strong ring exchange are respectively presented in Section 4.1 and 4.2. The nematic phase 
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in the regime of strong ferromagnetic bilinear coupling is discussed in Section 4.3. Section 5 
contains summary and concluding remarks. 

2. Model 



We consider the S = 1/2 Heisenberg spin ladder with four-spin ring-exchanges under a 
magnetic field. The model Hamiltonian is given by 

n = UQ + Uh (1) 

Hq = Jrung^ Si^i • S2,i + Jleg ^ (Si,r + S2,; • S2,i+l) 

+J4^(p, + pr') (2) 

I 

where s^^; is a spin-1/2 operator at a site on j-th leg and l-th rung and Jrung (<^icg) denotes 
the bilinear exchange constant on rungs (legs). The operator Pi {Pi~^) represents the four- 
spin ring-exchange process in which four spins on a plaquette {(1, 1), (2, 1), (2, 1 + 1), (1, 1 + 1)} 
exchange their positions clockwise (counterclockwise). Throughout this paper, we consider the 
case of Jrung = Jieg = J a^d J4 > and parameterize them as 

J = cos 6, J4 = sm9, (4) 

where < ^ < tt. 

In terms of the spin operators, the ring exchange is expressed in a complicated form. 

Pi + = Si,i • S2,l + Si,/+i • S2,/+l + Si,i • Si,i+i + S2,l ■ S2,/+l + Si,/ • S2,/+l + 82,; • Si,i+i 

+4(Si,/ • S2,/)(Sl,/+l • S2,Hl) + 4(Si,/ • Si,/+i)(S2,i • S2,i+l) 

-4(Si,/ • S2,i+l)(s2,i • Si,/+i), (5) 

where wc omit a constant term. Using this expression, the Hamiltonian (1) is rewritten in a 
form containing bilinear and biquadratic couplings. For the later use, we introduce an extended 
Hamiltonian with generalized four-spin exchanges defined as 

Hext = Si^l ■ S2,l + (Sl,/ • Si,/+i -|- S2,/ • S2,l+l) 

I I 
I 

-FJ„X (Si,/ • S2,/) (Si,/+i • S2,/+l) + ^11^(81,; • Si,;+i) {S2,l ■ 82,1+1) 
I I 

+Jdd^{Sl,l ■ S2,l+l) {S2,l ■ 
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-hY.Ki + 4,i)- (6) 

The ring-exchange model (1) with Eq. (4) is expressed as a parameter hne in the extended 
parameter space, Jj- = J + 2J4, Ji = J + J4, Jd = J4, and Jrr = Jn = —Jdd = 4J4. 

The ground-state properties of the ring-exchange model (1) have been studied for some 
limiting cases. For zero magnetic field h = 0, the ground-state phase diagram for the whole 
parameter region was determined by extensive numerical calculations.^^) It was shown that 
besides the conventional rung-singlet phase {6 < O.OTtt) and ferromagnetic phase (0.947r < 6) 
the model exhibited unconventional phases, including the staggered-dimer phase (O.OZtt < 6 < 
O.IStt), the scalar-chirality phase (O.IStt < 9 < O.SQtt), and the singlet phases with dominant 
vector-chirality (0.397r < 9 < O.SStt) and collinear-spin (O.SStt < 9 < 0.947r) short-range 
orders. All the phases except for the ferromagnetic one have a finite energy gap. Furthermore, 
the spin-chirality duality transformation,^^' which applies also to the case of /t > 0, revealed 
various duality relations between order parameters as well as duality mappings in the model 
Hamiltonian. We will discuss the duality transformation in the following section. 

For the antiferromagnetic two-leg spin ladders without ring exchanges, Eq. (1) with J > 
and J4 = 0, the effect of the magnetic field has been well understood.^^"^^) The energy gap 
above the singlet ground state dX h = decreases as the field h increases and closes at a 
critical field he- The system then enters a critical regime, where the low-energy excitations 
are described by the one-component Tomonaga-Luttinger (TL) liquid.^^'^^'^^^ Throughout the 
regime, the transverse spin correlation with {qx,Qy) = (tt, tt) dominates. The field-dependence 
of the TL-liquid parameter, which governs the low-energy physics of the system, as well as 
that of nonuniversal coefficients included in the TL-liquid theory were also obtained numeri- 
cally.^^' When h exceeds the saturation field hg, the fully polarized ground state emerges. 

Compared with the cases above, little is known about the case with J4 > and h > 0. 
Previous studies were performed for antiferromagnetic bilinear coupling J > and rather small 
ring exchange J4 < J. From numerical exact-diagonalization studies for small clusters with 
the level-spectroscopy analysis, it has been shown that for not too small J4 a plateau emerges 
in the magnetization curve at a half of the saturated magnetization.^^' The perturbation 
theory around the strong-rung-coupling limit Jrung — ^ 00 have concluded that the ground 
state at the plateau breaks the translational symmetry, resulting in the LRO of the staggered 
pattern of spin-singlet and spin-triplet rungs. A phenomenon called "77-inversion" , in which 
a longitudinal incommensurate spin correlation becomes stronger than the transverse-spin 
one, was predicted by exact diagonalization, while any direct observation of the behavior 
of the correlation functions has not been achieved yet.^^^ The possibility of the true LRO 
of vector-chirality was also discussed within the bosonization approach^°) though it is still 
elusive whether or not it realizes. 
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In this section, we discuss two duality transformations in the extended Hamiltonian Ticxt 
in Eq. (6) and the corresponding duahty relations between various phases. The one transfor- 
mation is the so-called spin-chirality duality and the other is its extension. The former was 
already introduced and discussed in Refs. 26 and 27 in detail, and we present here a brief 
summary for completeness. 

The spin-chirality duality transformation is defined as a mapping from original spin oper- 
ators in a rung, Si^i and S2,;, into new pseudospin ones; it is given by 

o 1 / 

= 2 l^l>' + ^2,/) - Sl,i X S2,l, 

S2,l = ^(si,i+S2,/)+Si,; XS2,/. (7) 

The new pseudospin operators obey the usual commutation relations for spins and satisfy 
= (^2,/)^ = 3/4. The spin-chirality transformation gives duality mappings between 
various order parameters. For example, the (g^:, 7r)-spin order operator, 

a(fe,7r) = ^e'^^'(si,,-S2,0, (8) 
I 

converts into the vector-chirality one, 

ac(g.) = 2j]e^''^'Si,^ xS2,/, (9) 

I 

and vise versa (up to an overall sign factor), while the total rung-spin operator Si_; -|- S2,; is 
self-dual under the transformation. We note that the duality transformation (7) corresponds 
to a unitary transformation,^^) 

Sn,i = u{{ei})sn,iu\{ei}), (10) 

with a unitary operator, 

U{m) = lle^P iei(^s,,i-S2,i-\'j , (11) 

and 01 = 7r/2 for arbitrary /. In this paper, we call this transformation the duality /. 

When applied to the extended model Ticxt, the duality transformation / leaves the form 
of the Hamiltonian unchanged and gives a duality mapping of the coupling parameters. To 
see the parameter mapping, it is convenient to rewrite the Hamiltonian as 

Wext = Si^l ■ S2,l + ^rr ^ (Sl,i • 82,;) (Si,;+i • S2,;+l) 

I I 

+W^{Si^l + S2,l) ■ (Si,,+i + S2,i+l) 
I 

+^ Yl [(^1'' ■ (^2,Z • S2,Z+l) + (Si,i • S2,i+l) (Si,i+i • S2,l)] 
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- S2,Z) • (Sl,Z+l - S2,i+l) + 4(si,i X S2,l) ■ (Si,i+i X S2,/+l)] 

I 

[(Si,i - S2,l) ■ (Sl,i+l - S2,i+l) - 4 (Si,i X S2,0 • (si,z+i X S2,;+l)] 

I 

+ (12) 

I 

where 

= ^(Jl + ^d), X=^iJn + JM), 

Z = -^(Ju- Jdd) + Jd). (13) 

After some algebra, one finds that the duahty transformation changes the sign of the cou- 
pling constant Z to — Z and leaves the other terms unchanged. Therefore, the duality / 
leads to the mapping in the seven-dimensional parameter space, ( Jj., Jj-r, VF, X, 1", Z, /i) to 
(Jr, Jrr, W, X, Y, —Z, h). In terms of the original coupling parameters, the mapping is given by 

J\ = 1^{J\ + Jd) + -^{Jl\- Jdd), Jd = -^{J\ + Jd) - -^{J\\- Jdd), 

Ju = 2(Ji- Jd) + ^(Jii + Jdd), Jdd = -2(Ji- Jd) + ^(^ii + Jdd)- (14) 

The extended model Hext with Z = 0, equivalently 4( Ji — J^) = Ju — Jdd, is therefore self-dual 
under the duality transformation /. 

Since the duality transformation / [Eq. (7) or Eq. (10)] consists of independent transforma- 
tions within each rung, one can construct arbitrary transformations with any set of the phase 
factors {Oi}. Among them, a set {9i = (— l)'7r/2} yields another useful duality transformation 
given by 

Sl,i = ^ (Si,, -I- S2,/) - (-l)'si,i X 82,^, 

S2,; = ^ {si,i + S2,/) + (-l)'si,, X S2,;. (15) 

We note that the transformation is a product of the duality / and the exchange of two spins in 
every second rung, Si / ^ S2^/ for, say, odd /. This transformation also yields another duality 
relations between order parameters; for instance, the (g^^, 7r)-spin order parameter OsigxiT^) 
is the dual of the vector-chirality one Ovc(fe + tt) with a momentum shift by vr. We call this 
transformation the duality //. 

Applying the duality transformation // to the extended model Ticxt causes a sign change 
on the y-term while the other terms are unchanged. The duality //therefore yields a mapping 
from (Jr, Jrr, W, X, Y, Z, h) to (Jr, Jrr, X, — Y, Z, h). The mapping of the original coupling 
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parameters is given by 



-'l = + Jd) - -^{Jn — Jdd), Jd = ^iJl^ Jd) + -g{J\\ - Jdd), 

Ju = -2(Ji- Jd) + ^(Jll + Jdd), ^dd = 2(Jl- Jd) + ^(^ll + ^dd)- (16) 

The extended model Hext in the parameter surface Y = 0, i.e., 4(Ji — Jj) = — (Ju — Jdd) is 
thus self-dual under the duality //. 

The duality mappings / and // enable us to make many exact statements on the phase 
diagram and the properties of the model. For example, if one knows properties of the model 
at a certain parameter point, one can immediately translate them to its three dual points. 
The mappings also yield a strong constraint on the phase diagram that the phase boundaries 
are symmetric with respect to the self-dual surfaces. Further, if there occurs a direct phase 
transition between two phases which are dual to each other, the transition line is exactly on 
the self-dual surfaces. Note that a self-dual phase, whose order parameter is self-dual under 
duality transformations / and/or II, can extend over the corresponding self-dual surfaces. 

When one considers the ring-exchange model (1), it should be noticed that neither the 
duality mapping / nor // preserves the form of the ring-exchange coupling. As a result, a 
parameter point on the ring-exchange model is mapped to a point which is not on the model. 
Even so, the duality relations give us useful results. The ring-exchange model (1) is self-dual 
at J = 2J4 [9 = Osdi = tan-i(l/2) ~ 0.1487r]26) and J = -2J4 [9 = 9^^!! = tan-i(-l/2) ~ 
0.8527r] under the duality mappings / and //, respectively. This is indeed consistent with 
the phase diagram aX h = obtained in Ref. 25 : The self-dual point I, J = 2J4, is the 
transition point between the staggered-dimer and scalar-chirality phases, which are dual to 
each other,^^) while a crossover between the singlet phases with the dominant staggered vector- 
chirality correlation and the (0,7r)-spin one occurs at the self-dual point //, J = — 2J4. Since 
the duality transformations do not change the magnetic field h, these dualities give two self- 
dual lines J = ib2J4 in the magnetic phase diagram in the 9 versus h plane. The presence 
of the self-dual line / and the fact that the model exhibits a TL liquid with a dominant 
(tt, 7r)-transverse-spin quasi-LRO for ^ = lead us to a prediction: If the region of the (tt, tt)- 
spin-dominant TL liquid extends up to the self-dual line /, there must be a crossover to its 
dual state, i.e., a TL liquid with the dominant staggered vector-chirality quasi-LRO, exactly 
at the self-dual line. We will see in the next section that this is indeed the case. Further, 
one may also expect that there occurs a crossover at the self-dual line // between the vector- 
chirality-dominant TL liquid for strong J4 and its dual state, (0, 7r)-spin-dominant TL liquid, 
for strong ferromagnetic J. As will be shown, this crossover sets in for a small field h while for 
a strong field a nematic phase, which is self-dual under the duality //, emerges and extends 
over the self-dual line //. 



7/18 




Fig. 1. (a) Magnetic phase diagram of the ring-exchange model (1) for < 6* < tt, where Jrung = 
Jieg = J = COS0 and J4 = sine. Symbols "P", "C^^", "VC^", "F" denote the regions of the 
magnetization plateau, transverse-spin-dominant TL liquid, vector-chirality-dominant TL liquid, 
and fully-polarized state, respectively. Circles represent the critical fields of the plateau bound- 
aries, hpi and /ip2, while the squares, diamonds, crosses respectively show the saturation field hs, 
boundaries he of the phases at zero magnetization, and the field /icusp at the cusp singularity. Solid 
curve show the saturation field obtained from the exact-diagonalization calculation (see text in 
Sec. 4.3). Vertical dashed lines are the self-dual lines / and //. The shaded area around 9 ^ 0.27r 
and small h indicates a region where the dominant correlation is not identified. Dotted curves 
are a guide for eye. (b) Enlarged phase diagram for 0.847r < 6 < 0.947r. Symbols "N" and "S" 
denote the nematic and spin-singiet phases, respectively, while triangles and diamonds show the 
boundaries of the phases. The bold crosses represent the boundaries of the region of two-magnon 
bound pairs, 6'mpi and 6'inp2- 



4. Numerical results 

In this section, we present our numerical results for the ring-exchange model Ti. [Eq. (1)] 
under a magnetic field. To determine the ground-state phase diagram, we have calculated the 
magnetization curve and various correlation functions. Our findings are summarized in Fig. 
1. The details of the results are discussed in the following. We first consider the region of 
< 0.847r, where we find the magnetization plateau and vector-chirality-dominant TL liquid. 
We then show our results on the nematic phase found for > 0.847r. 

4.1 Magnetization plateau 

We first discuss the magnetization curve. To obtain it, we have calculated the lowest energy 
Eq{M) of the Hamiltonian T^q [Eq. (2)] in a subspace characterized by the total magnetization 
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Fig. 2. Magnetization curve for (a) 9 = O.IOtt, (b) 6 = 0.307r, (c) 9 = O.SOtt, and (d) 9 = O.SOtt. The 
system size is Lj. = 51 and the number of the kept states in the DMRG calculation is up to 250. 
The horizontal dotted line represents the magnetization at the plateau M = (Lr + l)/(2Lr). Insets 
in each panel show the 1/Lr-dependence of the fields at the plateau boundaries, hpi and /ip2, and 
the saturated field hs- 



per rung, M = S^^^/Lj-, where Sfot = Z^/(sf 1 + ^21) -^r is the number of rungs. Note that 
M = 1 for the saturated state. We then determine the ground-state magnetization M{h) by 
searching for the magnetization M which minimizes the energy Eo{M;h) = Eq{M) — hS^^ 
for a given field h. The calculation was performed by the DMRG method^^"^^) for the finite 
systems with = 31, 41, and 51. (The reason to use the systems with odd Lr will be discussed 
later.) The open boundary condition was imposed for the sake of the efficiency of the DMRG 
method. 

Figure 2 shows the magnetization curves for typical values of 9 and Lj- = 51. It is 
clearly seen in the figure that the magnetization curve exhibits a plateau at the magneti- 
zation M = 1/2. [To be precise, for the open ladders with odd Lr the plateau appears at 
the magnetization M = (Lr -|- l)/(2Lr), which converges M — > 1/2 as Lr — > 00.] After the 
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Fig. 3. Local spin correlation at each rung, (si^j •S2,i), in the magnetization plateau at 6* = O.SOtt and 
Lr = 51. 

system-size extrapolation of the magnetic field at the upper and lower bounds of the plateau, 
hpi and hp2, we find that the plateau emerges for a wide range of 9, 6pi < 6 < 6p2, where the 
critical values are estimated to be Opi = (0.07ib0.01)7r and 0p2 = (0.83ib0.01)7r. The estimate 
of the lower critical value 6pi is consistent with the previous result ^pi = 0.069 obtained by 
the exact diagonalization.^^^ The extrapolated values of hpi and hp2 as well as the saturation 
field E^nd the energy gap at zero magnetization he are plotted in the phase diagram Fig. 1. 

It has been predicted by the perturbation theory around the limit of the strong rung- 
coupling^^) that the ground state in the plateau has the LRO of the staggered pattern of rung- 
singlet and triplet states. To confirm this, we have computed the ground-state expectation 
value of the correlation between the two spins in each rung, (si^/ -82,/). The results for 9 = O.bir 
are shown in Fig. 3. The data clearly demonstrate the appearance of the LRO in the plateau 
state: The rungs at the open edges are almost in the spin-triplet state and the staggered 
pattern of triplet and singlet rungs penetrates into the bulk without a decay. Therefore, the 
translational symmetry is broken spontaneously in the plateau and the system has the two-fold 
degenerate ground states, being consistent with the necessary condition for the presence of the 
plateau.^^) We note that one of the two ground states is selected in the DMRG calculation due 
to the open boundaries. (This is the reason to select odd in the analysis of magnetization 
plateau; for even Lj-, there appears a kink between two different patterns from the open 
boundaries at the center of the ladder, resulting in an artificial step in the plateau of the 
magnetization curve.) The appearance of the LRO is observed for the whole parameter regime 
of the magnetization plateau. 

Another feature to be noted for the magnetization curve is a cusp singularity appearing 
at a magnetization M < 1/2. [See Fig. 2 (a) and (b).] The cusp is found for 6pi < 9 < 
0.57r, which corresponds to the region where the plateau exists and the bilinear exchanges 
are antiferromagnetic. We note that such a cusp singularity, which suggests a change of the 
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Fig. 4. Absolute values of correlation functions for (a) {9,M) = (O.Itt, 0.4), (b) (O.Itt, 0.7), (c) 
(O.Stt, 0.4), and (d) (O.Stt, 0.7). Squares, triangles, and circles represent the correlation functions 
C^^{r), C^^{r), and C^Q{r), respectively. The system size is — 60 and the number of the kept 
states in the DMRG calculation is up to 300. Truncation errors are smaller than the size of the 
symbols. 

structure of the excitation spectrum, is often observed in frustrated spin systems. ^'^^ The 
strength of the field at the cusp, /icusp, is also plotted in the phase diagram, Fig. 1. 

4-2 Vector chirality 

Besides the gapped phases of the zero magnetization and the magnetization plateau, there 
is a wide region of critical states in the phase diagram. To characterize the critical states and 
identify the dominant order, we have calculated various two-point correlation functions in the 
ground state of the finite open ladders with up to Lj- = 60 rungs using the DMRG method. The 
correlation functions considered are the ones of total {qy = 0) rung spin s"i + {a = x,z), 
Neel-type {qy = vr) rung spin s" ^ — S21, vector chirality (su x S2,;)°, scalar chirality Si_; • (s2,i x 
Si^/_i_i), staggered dimers Si^; • Si^i^i — S2^i ■ S2,i+i, and two-magnon pairing operators 
and s^iS^i^i- To lessen the open boundary effects, we calculated the correlation functions on 
six different pairs of two rungs (l, I') for each distance r = | / — /'|, selecting the center position 



11/18 



J. Phys. Soc. Jpn. 



Full Paper 



.(a) 



- ■ ■ ■ 


A 


A 


A 


A 


A 


A i 


- ■ ■ ■ 














- ■ ■ ■ 


A 


A 


A 


A 


A 


A i 


- ■ ■ ■ 


A 


A 


A 


A 


A 


A i 


■ ■ ■ o 


o 


o 


o 


o 


o 


o c 


- ■ ■ ■ 




A 


A 


A 


A 


A i 


- ■ ■ ■ 




A 


A 


A 


A 


A i 


■ ■ X 




A 


A 


A 


A 


A d 


' ■ ■ X 






A 


A 


A 


A i 



A 


A 




♦ 


♦ 


A 


A 




♦ 


♦ 


A 


A 


« 


♦ 


♦ 


A 


A 




♦ 


♦ 


o 


A 




♦ 


♦ 


A 


A 




♦ 


■ 


A 


A 




♦ 


■ 


A 


A 




■ 


■ 


A 


A 




■ 


■ 



Fig. 5. Dominant correlation functions for (a) < < O.Stt and (b) O.Stt < 9 < O.Ott. Squares, tri- 
angles, and diamonds represent the critical states where the transverse Neel-type-spin, transverse 
vector-chirality, and magnon-pairing correlation functions are dominant, respectively. Open circles 
show the plateau state. Vertical dashed lines are the self-dual lines / and //. Crosses denote the 
point where boundary effects are too strong to determine the dominant order. 



^0 = + ^0/2 to be as close as possible to the center of the ladder . We then took the average 
as an estimate of the correlation function for each distance r. 

Among the correlation functions considered, the transverse Neel-type-rung-spin and 
vector-chirality correlation functions, 

C!^{r) = j^{{sli-sli){sl,-sl,)), (17) 

C^cir) = ^((si,,xs2,,r(si,,'XS2,rr), (18) 

turn out to be dominant in the critical phase for 9 < 0.847r. In Fig. 4, we show the data of 
these correlation functions with the longitudinal total-rung-spin one, 

C!oir) = ^ [{{sli + sl,){sl, + 4,)) - {4,1 + sli){sl, + si,)] . (19) 

The Neel-type rung-spin correlation C^(r) and vector-chirality one C^^{r) decay algebraically, 
oscillating with the wavenumber = tt. Since the decay exponents of these correlations seem 
to be the same, we determine which is the dominant one by comparing their amplitudes. 
Figure 5 summarizes the results in the 9 versus M plane. The Neel-type rung-spin correlation 
C^(r) is dominant for 9 < ^sdi while the vector-chirality one C^^{r) is for 9 > ^sdi^ The 
crossover occurs at the self-dual line 7, as discussed in Sec. 3. The difference between the 
amplitudes of C^(r) and C^^{r) becomes smaller as the magnetization M increases and they 
converge the same value at M — 1 (see Fig. 4). The vector-chirality-dominant TL liquid thus 
emerges in the wide parameter space of strong ring exchange including both the regions of the 
antiferromagnetic and ferromagnetic bilinear exchanges. We note that there remains a small 
region around 9 ~ 0.2tt and small M where too strong open boundary effects prevent us from 
determining the dominant correlation functions. For clarifying the nature of the model in the 
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Fig. 6. (a) Magnetization curve for 9 = O.SStt. The system size is = 51. (b) Field dependence of 
excitation energy Eq{M; h) — £'gs(/i). Solid and dotted lines represent the energy level with even 
and odd S^^^, respectively. 



region, it is required to perform the calculation for much larger systems, which is out of the 
scope of our numerics and left for future studies. 

We mention the possible "r/-inversion" in the system. For the present model (1), it was 
predicted that around the plateau there is a parameter regime where the longitudinal in- 
commensurate spin correlation function becomes stronger than the transverse staggered spin 
Qj^g 39) have found that for M > 0.5 the decay exponent of the longitudinal total-rung- 
spin correlation function CgQ(r) approaches those of Cf^(r) and C^f.{r) as M decreases and 
they seem very close to each other at M = 0.6. This seems consistent with the prediction. 
However, unfortunately, we do not obtain a clear evidence of the phenomenon because of 
the incommensurate character of CgQ(r) as well as the open boundary effects, which make it 
difficult to estimate the decay exponent of C^q precisely. To overcome the difficulties, numer- 
ical calculations for larger systems and, more preferably, theoretical descriptions such as the 
bosonization analysis for the low-energy states of the system would be necessary. 

4-3 Nematic phase 

Here, we discuss our results for 6 > 0.847r, where we find the nematic phase. Figure 6 
(a) shows the magnetization curve obtained by the DMRG calculation for 6 = O.SSvr. An 
interesting feature is that for a field h larger than a certain critical value the magnetization 
changes by a step A5^^q^- = 2. This suggests that magnons in the background of the saturated 
state form bound pairs of two magnons. To confirm this, we calculate the field dependence of 
the excitation spectrum of the lowest level in each subspace of M, Eq{M; h) — Egs{h), where 
Egs{h) is the ground-state energy under a field h. The result in Fig. 6 (b) clearly indicates 
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Fig. 7. Absolute values of correlation functions for {9, Ad) = (O.SStt, 0.7). Diamonds, circles, squares, 
and triangles represent the correlation functions Cyy^{r), CgQ(r), C^^(r), and C^^{r), respectively. 
The system size is L,. = 60 and the number of the kept states in the DMRG calculation is up to 
300. Truncation errors are smaller than the size of the symbols. 

that for large h only the levels with odd S^^^ becomes the ground state while the levels with 
even S^^^ always have a finite excitation energy. We thereby conclude that the magnons indeed 
form two-magnon bound pairs with a finite binding energy. [Note that, since is odd, even 
(odd) magnons are included in the state with odd (even) Sfo^.] 

Since the magnetization curve in Fig. 6 was calculated for the system with the open 
boundary condition, one may think that the step by A54t = 2 is due to open boundary 
effects. To answer the question and determine the boundary of the region of the magnon pairs 
more accurately, we have performed the exact diagonalization calculation for the periodic 
ladder within the subspace of n-magnon states, where the number of magnons is up to n = 4. 
The system size treated was up to = 48, which turned out to be large enough. We calculate 
the lowest energy in each subspace and analyze the magnon instability in the saturated state. 
We thereby find that for 9 < ^mpi = 0.8447r the single-magnon instability is the strongest and 
the saturation field is given by hs = 3J + 8J4 = 3 cos + 8 sin ^, which is obtained analytically. 
However, when 9 exceeds ^mpij the two-magnon instability dominates and determines the 
saturation field h^. The region of two-magnon bound pairs extends for ^mpi < 9 < 0mp2 = 
0.9237r. In the small region between 0mp2 and the boundary of the ferromagnetic phase 9 = 
0.9357r, we find that four-magnon instability takes place, but we must note that the instability 
of more-magnon bound pairs may emerge in the vicinity of the ferromagnetic phase. The 
saturation field obtained and the boundaries of the two-magnon pairing region, ^mpi and 
0mp2) are plotted in Fig. 1. 

To investigate the nature of the system in the magnon-pairing regime in more detail, we 
have calculated the several correlation functions mentioned above using the DMRG method. 
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From the calculation, it turns out that the magnon-pairing correlation functions, 




(20) 



(21) 





(22) 



dominate in the region where the two-magnon bound pairs are found. A typical example 
of the data of Cmp(r) are shown in Fig. 7 with the longitudinal total-rung-spin correlation 
CgQ{r), the transverse Neel-type rung-spin one Cg^{r), and the transverse vector-chirality 
one Cy(,(r). The magnon-pairing correlation function C{S'p{r) is clearly the strongest: CgQ(r) 
decays algebraically but faster than Cmp(r) while C^^{r) and C^^{r) decay exponentially. 
We note that the other magnon-pairing correlation functions C^p(r) and Cmp(r) also decay 
algebraically with the same exponent as that of Cmp(r) while the amplitudes are different. 
Interestingly, the magnon-pairing correlation functions have the sign Cmp(r) > 0, C^^p(r) > 
0, and Cmp(r) < for arbitrary distance r; they exhibit the "d-wave-like" symmetry. All 
these behaviors are in common with those of the nematic phase found for the ring-exchange 
model in the square lattice,^^^ in which phase the two-magnon bound pairs undergo a bose 
condensation resulting in the magnon-pairing order, though in our one-dimensional system the 
true LRO of the pairing correlation is destroyed by a strong quantum fluctuation. We therefore 
conclude the ring-exchange ladder (1) for the parameter regime is in (one-dimensional analog 
of) the nematic phase. The region of the nematic phase, where we find the magnon-bound 
pairs and the dominant magnon-pairing quasi-LRO, are shown in Figs. 1 (b) and 5 (b). The 
phase appears at the fields h larger than a critical value hmp and extends over the self- 
dual line //. This is allowed since the magnon-pairing order parameter is self-dual under 
the duality transformation. The decay exponents of the magnon-pairing correlation Cmp(r) 
and longitudinal total-rung-spin one CgQ(r) become closer to each other as the magnetization 
decreases and seem almost the same for ^ = 0.85 and M = 0.3. 

At a field h lower than the critical value /imp , we find that the vector-chirality correlation 
Cy^,(r) is dominant for 9 < ^sdii while the Neel-type-rung-spin one C^^(r) with wavenumber 
qx = along the leg direction dominates for 6 > Ogdii- The crossover between the vector- 
chirality-dominant and Neel-type-spin-dominant TL liquids occurs at the self-dual line //, as 
expected in Sec. 3. 

5. Summary and Concluding remarks 

In summary, we have studied the spin- 1/2 two-leg Heisenberg ladder with four-spin ring 
exchanges under a magnetic field, for a wide parameter regime including both the antifer- 
romagnetic and ferromagnetic bilinear couplings. We have introduced a duality transforma- 
tion, which is a simple extension of the spin-chirality duality transformation developed previ- 
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ously^^'^^^ but leads to a nontrivial duality mappings on coupling parameters in the Hamil- 
tonian as well as order parameters. These dualities yield two self-dual lines in the parameter 
space of the ring-exchange ladder (1). 

Further, using the density-matrix renormalization-group and exact diagonalization meth- 
ods, we have determined numerically the magnetic phase diagram. Fig. 1, including the magne- 
tization plateau at Af = 1/2 and the regions of TL liquids with different dominant quasi-LRO. 
We have found the vector-chirality-dominant TL liquid emerging in a wide parameter regime 
of the strong ring-exchange coupling between the self-dual lines, while there appear TL liquids 
with the dominant transverse spin quasi-LRO outside of the self-dual lines. Moreover, we have 
identified the nematic phase, in which the magnons form bound pairs and the condensation 
of the bound pairs leads to the dominant magnon-pairing correlation functions, in a finite 
regime in the vicinity of the ferromagnetic phase. 

The formation of the magnon bound pairs and their condensation in the nematic phase 
are not peculiar to the ring-exchange ladder (1) but are found in various S = 1/2 systems: 
The phenomena have been reported for the ring-exchange model in the square lattice'^^-' and 
triangular lattice, and the zigzag ladder with ferromagnetic nearest-neighbor and antifer- 
romagnetic next-nearest-neighbor couplings. ^^"''^'^^ These observations suggest that the phe- 
nomena are a general feature of "frustrated ferromagnets" in which the fully-polarized state 
is destabilized by certain perturbations such as ring exchange and/or antiferromagnetic cou- 
plings. Systematic studies to develop theoretical descriptions of the state and to clarify the 
relation between the nematic phases in 5' = 1/2 and S* > 1 systems would be important for 
understanding the phenomena. 

For the two- leg ladder with extended four-spin exchange Hcxt [Eq. (6)], it was predicted 
in Ref. 40 that the model under a magnetic field ft, > could realize the true LRO of the 
longitudinal vector-chirality (si ^ x S2^;)^, while the author also pointed out that the form of 
the ring-exchange coupling, J„ = J\\ = — Jdd) was not suitable for the appearance of the LRO. 
Indeed, we have not observed the appearance of the true LRO of the vector-chirality in our 
calculation on the ring-exchange ladder (1). Searching for the LRO in a ladder model with 
the generalized four-spin exchanges of more suitable form would be an interesting issue. 
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